Vector Meson production from NLL BFKL 
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The amplitude for the forward electroproduction of two light vector mesons can be 
written completely within perturbative QCD in the Regge limit with next-to-leading 
accuracy, thus providing the first example of a physical application of the BFKL ap- 
proach at the next-to-leading order. We study in the case of equal photon virtualities 
the main systematic effects, by considering a different representation of the amplitude 
and different optimization methods of the perturbative series. 

1 Introduction 

In the BFKL approach pQ, both in the leading logarithmic approximation (LLA), which 
means rcsummation of all terms (a s ln(s)) n , and in the next-to-leading approximation 
(NLA), which means resummation of all terms a s (a s ln(s))™, the amplitude for a large- 
s hard collision process can be written as the convolution of the Green's function of two 
interacting Reggeized gluons with the impact factors of the colliding particles. 

The Green's function is determined through the BFKL equation. The kernel of the BFKL 
equation is known now both in the forward [2] and in the non-forward [3] cases. On the 
other side, impact factors are known with NLA accuracy in a few cases: colliding partons [2], 
forward jet production [5] and forward transition from a virtual photon 7* to a light neutral 
vector meson V = p°,uj,<fi [BJ. The most important impact factor for phenomenology, the 
7* — > 7* impact factor, is calling for a rather long calculation, which seems to be close to 
completion now pj [8] . 

The 7* — > V forward impact factor can be used together with the NLA BFKL forward 
Green's function to build, completely within perturbative QCD and with NLA accuracy, the 
amplitude of the 7*7* — ► VV reaction. This amplitude provides us with an ideal theoretical 
laboratory for the investigation of several open questions in the BFKL approach. Besides, 
this process can be studied experimentally at the future at ILC, see Refs. [5J. 

2 Representations of the NLA amplitude 

The process under consideration is the production of two light vector mesons (V — p , w, </>) 
in the collision of two virtual photons, 7*(p) 7*(p') — * V(pi) Vfa). Here, neglecting the 
meson mass my, p\ and P2 are taken as Sudakov vectors satisfying pi = p\ = and 
2(pif>2) = s; the virtual photon momenta are instead p = api — Q\/ (as)p2 and p' = 
a'p2 — Q\l '(ct' s)pi, so that the photon virtualities turn to be p 2 = —Q\ and (p') 2 = — Q%. 
We consider the kinematics when s 3> Q\ 2 3> ^qcd an< ^ a = 1 + Q\l ' s + C( s_2 ): a' = 1 + 
Q 2 /s + 0(s~ 2 ). In this case vector mesons are produced by longitudinally polarized photons 
in the longitudinally polarized state [BJ. Other helicity amplitudes are power suppressed, 
with a suppression factor ~ my/Qi.2- We will discuss here the amplitude of the forward 
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scattering, i.e. when the transverse momenta of produced V mesons are zero or when the 
variable t = (pi — p) 2 takes its maximal value to = — Q\Q\j s + 0(s~ 2 ). 

The NLA forward amplitude can be written as a spectral decomposition on the basis of 
eigenfunctions of the LLA BFKL kernel: 
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Here the bulk of NLA kernel corrections are exponentiated, a s = a s N c /ir and D12 = 
— 4:ire q fv/(N c Qi,2), where fv is the meson dimensional coupling constant (f p w 200 McV) 
and e q should be replaced by e/y/2, e/(3-\/2) and — e/3 for the case of p°, oj and (j> meson 
production, respectively. Two scales enter the expression ([1]), the renormalization scale 
and the scale for energy s - 

Alternatively, the amplitude can be expressed as a series: 

Q1Q2 Im s (Aseries) 1 / N 2 /n\ 
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The b n coefficients are determined by the kernel and the impact factors in LLA, while the 
d n coefficients depend also on the NLA corrections to the kernel and to the impact factors. 
We refer to Ref. [TD] for the details of the derivation and for the definition of the functions 
entering these expressions. 



3 Numerical results 

In Ref. [TU] we presented some numerical results for the amplitude given in Eq. Q for 
the Qi = Q2 = Q kinematics, i.e. in the "pure" BFKL regime. We found that the d n 
coefficients are negative and increasingly large in absolute values as the perturbative order 
increases, making evident the need of an optimization of the perturbative series. We adopted 
the principle of minimal sensitivity (PMS) by requiring the minimal sensitivity of the 
predictions to the change of both the renormalization and the energy scales, fiR and Sq. 
We considered the amplitude for Q 2 =24 GeV 2 and n/ = 5 and studied its sensitivity to 
variation of the parameters fj,R and Yq = ln(so/Q 2 )- We could see that for each value of 
Y = ln(s/Q 2 ) there are quite large regions in [ir and Yq where the amplitude is practically 
independent on [ir and Yq and we got for the amplitude a smooth behaviour in Y (see 
the curve labeled "series - PMS" in Figs. [1] and [21). The optimal values turned out to be 
/i/f ~ 10Q and Yq ~ 2, quite far from the kincmatical values fin = Q and Yq — 0. These 
"unnatural" values probably mimic large unknown NNLA corrections. 



DIS2007 




As an estimation of the systematic effects in our determination, we considered also 
the "exponentiated" representation of the amplitude, Eq. (TTJ), and different optimization 
methods. For more details on the following, see Ref. [12] . 

At first, we compare the series and the "exponentiated" determinations using in both 
case the PMS method. The optimal values of /j,r and Yq for the "exponentiated" amplitude 
are quite similar to those obtained in the case of the series representation, with only a slight 
decrease of the optimal Fig. Q] (left) shows that the two determinations are in good 
agreement at the lower energies, but deviate increasingly for large values of Y. It should 
be stressed, however, that the applicability domain of the BFKL approach is determined by 
the condition q. s {)jlb)Y ~ 1 and, for Q 2 =24 GeV 2 and for the typical optimal values of /in, 
one gets from this condition Y ~ 5. Around this value the discrepancy between the two 
determinations is within a few percent. 

As a second check, we changed the optimization method and applied it both to the series 
and to the "exponentiated" representation. The method considered is the fast apparent 
convergence (FAC) method [13] , whose strategy, when applied to a usual perturbative ex- 
pansion, is to fix the renormalization scale to the value for which the highest order correction 
term is exactly zero. In our case, the application of the FAC method requires an adaptation, 
for two reasons: the first is that we have two energy parameters in the game, /ir and Yq, the 
second is that, if only strict NLA corrections are taken, the amplitude does not depend at 
all on these parameters. For details about the application of this method, we refer to [12]. 
Here, we merely show the results: the FAC method applied to the series representation (see 
Fig. [T] (right)) and to the exponentiated representation (see Fig. [2] (left)) gives results in nice 
agreement with those from the PMS method applied to the series representation, over the 
whole energy range considered. 

Another popular optimization method is the Brodsky-Lepage-Mackenzie (BLM) one [H] , 
which amounts to perform a finite renormalization to a physical scheme and then to choose 
the renormalization scale in order to remove the /3o-dependent part. We applied this method 
only to the series representation, Eq. ([2J). The result is compared with the PMS method in 
Fig. [2J (right) (for details, see Ref. [12] V 
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Figure 2: Im 3 {A)Q 2 /(s DiD 2 ) as a function of Y at Q 2 =24 GeV 2 (n f = 5): (left) series represen- 
tation with PMS and "exponentiated" representation with FAC, (right) series representation with 
PMS and with BLM. 
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